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ABSTRACT 


An analysis of test system misalignment is presented for dynamic tension 
test. Sheet type rectangular 1100-0 aluminum specimens are used for discussion. 

For a constant strain rate tension test, the strain rate is constant only 
on the neutral axis of the specimen. The following results have been obtained: 
(a) The lower the strain rate is, the more significant the misalignment errors 
become, (b) Misalignment errors of 50% have been found at the extreme fibers 
of the specimen, (c) The strain rate variation in the cross-section decreases 
with increasing plastic strain and vanishes at plastic strain equal to 0.8% at 
the midspan of the specimen, (d) The neutral axis will shift toward the 
centerline of the specimen as the plastic strain increases. But, it will 
reach a limit and will not completely move back to the centerline. 

A more restricted uniform strain rate formulation is also presented. 

The result is compared with that of the nonuniform strain rate formulation. 


] . INTRODUCTION 


The effect of misalignment In a tensile test system has been previously 
investigated by a number of Investigators [1,2]. These Investigators have 
also Identified the sources of misalignment as poor conformance of specimen 
centerline to top and bottom grip centerlines, poor alignment of the top and 
bottom grip centerlines, and inaccurate machining of the test specimen Itself. 

A more detailed description of test system misalignment and the discussion 
concerning the need for such investigation may be found in Refs. [1,2] and also 
in Wu and Rummler [ 3 ] . 

In [3], a comprehensive analysis of misalignment effect was presented 
for tension test under static loading condition. In the analysis, three cases 
of misalignment were considered, which were the symmetric case, the cantilever 
case, and the case of the pinned specimen with eccentricity. It was found 
that the symmetric case is the most critical one as far as the effect of mis- 
alignment is concerned. In addition, Wu and Rummler concluded that the stress- 
strain curve is significantly affected by misalignment at strain levels corre- 
sponding to the knee portion of the stress-strain curve. Moreover, the strain 
at the outermost fibers of the specimen is strongly affected by misalignment 
and the misalignment effects are smaller for load trains with longer pull rods. 

The dynamic mechanical properties of materials have become increasingly 
important in the recent years. Tension tests have been conducted at various 
fixed strain rates to determine the dynamic behavior of materials. The problem 
of test system alignment which exists in the static test would also carry over 
to the dynamic test. However, to the present investigators' knowledge, the 
alignment problem under dynamic tensile load has not been investigated in the 


literature. 
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In this paper, an analysis of test system misalignment is presented for 
dynamic tension test. For simplicity, rectangular specimens of the sheet type 
ore considered in this investigation. Only the symmetric case of misalignment 
is studied due to the flnd.tng of [3] that this is the most severe case. Speci- 
mens of 1100-0 aluminum have been chosen for investigation. The reason for 
using this material as an example for discussion is that the constitutive 
equation employed in the present analysis has been shown to describe both the 
dynamic stress-strain behavior and creep for this material [4] and the material 
constants have already been determined. Theoretically speaking, the method 
presented herein would also apply to other materials. 

The purpose of this study is to investigate the effect of strain rate on 
the misalignment problem. The equations derived are used to calculate the 
strain rate distribution and the misalignment error at the center section of 
the specimen. The center section is of Interest, since the strain measuring 
devices (extensometer , strain gage, etc.) are usually attached to the center 
of the specimen. 

In Section 2, the constitutive equation used in the analysis is briefly 
described. The field equations for the misalignment problem found in the 
dynamic tension test and their discretized numerical models are formulated in 
Section 3. In the formulation, the strain rates are nonuniform across a 
cross-section of the specimen although the specimen is subjected to a constant 
strain rate test (constant strain rate at the neutral axis (N.A.)). A more 
restricted formulation where strain rate is constant all over the specimen is 
included in Appendix A. The results of this uniform strain rate formulation 

& 

The ultimate goal of this research progam is to determine the effect of mis- 
alignment in the creep test. The results obtained in this report will be 
used directly in the creep investigation. 


3 


serve as the Initial guess for obtaining numerical solutions of the highly 
nonlinear equations of Section 3. A detailed computing procedure for solution 
of this set of nonlinear equations Is described in Appendix B, 

The numerical results and discussion are presented In Section 4. It is 
believed that the effort made in this research contributes toward a more 
thorough understanding of the misalignment effect In the tension test. 
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2. THE CONSITUTIVB EQUATION 

The constitutive equation used in this investigation is the same as that 
given by Wu and Yip [5] for strain rate sensitive materials. This is a modified 
version of the endoehronlc constitutive equation originally proposed by Volanis 
[6]. For constant strain rate monotonic tension test, the constitutive equation 
is 

cr - ~ (1 + 65) [a Q - (o Q - o y ) (1 + 6C)“ n ) (1) 


in which a is the stress; o Q is the intercept of the asymptotic line of the 
reference stress-strain curve with the stress axis; o y is the yield stress; 

6 and n are material constants; z is the Intrinsic time defined by 

5 * k(0)0 (2) 


and k is the strain rate sensitivity function governed by 


k(6) 



(3) 


where 

O = e - o/E ( 4 ) 

is the plastic strain, e is the total strain, E is Young's modulus, 6 is a 

s 

material constant, and 0^ is a reference strain rate. 

Equation (1) represents a set of constant 0 stress-strain curves. For 

1100-0 aluminum at 150°C, the function k(0) was determined by Wu and Chen [4] 

-4 3 -1 

from a set of experimental data for the strain rate range of 10 to 10 sec 
This function is reproduced in Fig. 1, where 6 and 0 are equal to 

S <L\ 


-2 -5 -1 

3.643 * 10 and 1.30 x 10 sec , respectively. Figure 2 shows three a vs. 

0 curves, each corresponding to a constant 0 with corresponding k values equal 
to 1,00, 0.75 and 0.60. These three curves are used in the study of misalign- 
ment to be reported in the next section. 

Finally, it must be mentioned that equation (1) describes only the plastic 
behavior of the material under consideration. 
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3. FORMULATION OF THE MISALIGNMENT PROBLEM 

When a misaligned specimen undergoes a constant strain rate dynamic 
tensile test, it is expected that the strain rate varies both along the 
lateral and the longitudinal directions of the specimen. The strain rate 
remains constant only along the N.A. which is defined as the location of the 
mean stress (load/area) os in Ref. [3], 

Specimen dimensions and the coordinate system ore chosen the same way 
as in Ref. [3], and the latter Is shown in Fig. 3(a), Figure 3(b) shows the 
load train configuration of the symmetric case of misalignment. In the 
analysis, it Is assumed that plane cross-sections remain plane during defor- 
mation. Thus, the strain e at a point in the specimen is given by 

*<*’*> * 5 - 7m <5> 


where x and y are coordinates shown In Fig, 3(a), e is the strain on N.A., 
and p(x) is the radius of curvature of the specimen. The radius of curvature 
Is expressed in terms of the strains of the top fiber and the bottom fiber 
Eg through the expression 


p(z) 



( 6 ) 


in which h Is the width of the specimen. The distance d(x), which specifies 
the location of the N.A. and is measured from the bottom fiber of the specimen 
is given by 

d(x) e B " E 

=3 — — 

h e„ - e, 


( 7 ) 
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Due to Che variation of strain rate In the specimen, the misalignment 
problem is quite complicated mathematically* Only numerical solutions are 
feasible, Therefore, the specimen la discretized as shown In Fig. 4 . The 
stress, the strain, and the strain rate in each element are assumed to have 
constant values. Only half of the specimen needs to be considered due to 
the symmetric loading condition. Thus, the half span is cut into I sections 
along the longitudinal direction and resulted in I + 1 planes to be investi- 
gated (see Fig. 4 a). The lateral direction (the width h) is cut into J 
sections as shown in Fig, 4 b. 

Based on the discrete model and the constitutive equation ( 1 ), the 
stress-strain relation in the element (i,j) is 

0 ij " (1 + 05 ij^ [ff 0 " <*0 ~ V^ 1 + 8c ij^ 5 

with 1 ** 1, . . . ,1+ 1 5 j * 1,.,,,J (8) 

While equation ( 1 ) is exact, an additional assumption is made to arrive at 
equation (8). It is assumed that a generic point in the specimen maintains 
the same strain rate throughout the test although the strain rate may vary 
from point to point in a cross-section. Since the strain rate at a generic 
point varies only very slightly (this is justified by the numerical results), 
this assumption constitutes a close approximation when equation (8) is used 
to compute the stress. 

The strain rate sensitivity function is 



( 9 ) 
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Since the reference plastic strain rate d R may be different than Che strain 
*!* 

rate at the N»A», u, it is convenient to rewrite equation ( l J) as 


k lj ■ 1 - “a ln 


ff) ■ e « *" (t) 


( 10 ) 


It is remarked that 0 is prescribed for each test and 0 R is known from tile 
experimental data (see Ref. [4]). Hence the last term on the right-hand side 
of equation (10) is determined. The strain race sensitivity function on the 
N.A., denoted by k, may be obtained by putting 0^ * 0 in equation (10). 

Thus , 


k * 1 - 8 £n 
s 


(i\ 

VJ 


(id 


and equation (10) reduces to 


- k - 0„ £n 

0 / 


ij 




(12) 


In the numerical computation, the strain rates are expressed by the increments, 
i.e. , 


k y - k - S s in 


&) 


(13) 


Furthermore, from equations (4) and (2), it is obtained 


■ c « ■ °« /E 


(14) 


and 


- k u (E ij - V® 


(15) 
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The linear distribution of strain in ? ' , ro 88 ~section renders 


'ij IT 


+ " e iT ) (“ J T ■'-) ! 


.} ** 1| Mt |J 


(16) 


for the strain at the 1th plane (see Pig. 4b), The strains e iT 
referred to the continuous, deformed cross-section, whereas 
to the discretized planes of the numerical model. 

Balance of force acting on the specimen gives 


and e^jj are 
is assigned 


P 


er dA * b 

k 


a dy 


(17) 


In which P is the force, b is the thickness and A Is the cross-sectional area 
of the specimen. In a discrete form, this equation leads to 


bh 

J 


U 


ij 




for the ith plane. The expression h/J is the width of the element under con- 
sideration at the ith section. 

The stress a and the Intrinsic time 5 at the N.A. are related to the 
load P by 


j - a - ~ (1 + Be) [ffQ- (Oq “ o y )(l + 3i)~ n ] 


(19) 


Note that the quantities C, k, and a are assigned as input during numerical 
calculation. 

The balance of moment is given by 


M(x) 


Ja 


a y dA 


psP + bp 


at dy 


('V'N 
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In tha discrete model, this equation bccomoB 


M. 


Pi e P + 


bh 

J 


J 

I 

J*1 




t * • * 


, 1+1 


( 21 ) 


for the 1th plane. 

Referring to Fig. 3(b), the following equation is obtained 


M 


1 


- P 



+ u. 



( 22 ) 


where 

6* - 6 + (d - h/2) (23) 

and Is the deflection at ith section, 6 is the eccentricity, L is the pull 

rod length, and d denotes the location of the N.A, at the grip end. The 
s 

distance d„ is that of d evaluated at 1 * I + 1 and is further given by 
s 


^14*1 B ~ ^ 

d » h 

8 e I+l,B ~ ®I+1,T 


(24) 


where ^ and are strains of top and bottom fibers respectively at 

the grip end. 

The curvature relation is 


.C%, 

;V 


JL 



(25) 


which in its difference form is 


E iB g iT 
h 


u i+l 




i - 1,.. .,1- 1 


(26) 
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where Ax is the uniform distance between cross-sectional planes. 

The formulation of the misalignment problem under dynamic loading 
condition Is now complete. The set of equations (8), (15), (18), (21), (24) and 
(26) constitutes a system of highly nonlinear equations, which can be solved 
numerically for variables o^, e iT> e iB » 5^ » u i and d g for a specified load 
P subjected to the following boundary conditions: 


du 

dx 


at 


0 


(27) 


and 

u “ 0 at x = A/2 (28) 

In the above equation, A is the length of the specimen. 
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4. RESULTS AND DISCUSSION 

During analyses, three different; values of strain rate are assigned to 

the N.A. ns Input data to Investigate the Influence of strain rate, on misaligned 

“5 -2 

tensile specimen. The strain rates investigated are 1.30 * 10 , 1.24 x 10 , 

-1 -1 

and 7.63 * 10 sec and their corresponding values of strain rate sensitivity 
function k are 1.00, 0 . 7 !> and 0.60. 

The sheet type rectangular speclments are used. The geometrical data for 
the specimen are 

pull rod length/specimen length ■ 2 
specimen length (Jl) * 63.5 mm 
specimen width (h) = 12.7 mm 
specimen thickness (b) = 3.2 mm 
eccentricity 6 « 0.05 h 

In order to reduce computer time, only four sections have been cut along 
both the longitudinal and the lateral directions of the specimen, i.e., 1=4 
and J = 4. Thus, the unknown variables have been reduced to 55 In number and 
the set of nonlinear equations discussed In Section 3 consists now of 55 
equations. Due to the minimum number of elements in the discrete numerical 
model, some unstable numberical results are experienced at the grip end at 
the larger strain level and will be discussed later In the text. Nevertheless, 
this model provides a consistent output under a variety of loading conditions. 
Since the area of interest for the misalignment analysis lies in the midspan 
of the specimen (the measuring devices are usually attached to the midspan) , 
most of the results to be discussed are related to this area. 
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The normalized strain rate sensitivity function k/k for the top and 
bottom fibers at the midspan Is plotted In Fig. 5 against the plastic strain 
0 for the three strain rate levels previously mentioned. It Is seen that k 
decreases (increases) rapidly at small strain level for top (bottom) fiber 
until 0 reaches 0.8%, when k becomes uniform across the cross-section. Hence, 
the variation of strain rate In a cross-section diminishes as strain Increases. 
This phenomenon can be explained from the following equation 



which is obtained by rearranging equation (13). At the small plastic strain 
range, k/k is greater (smaller) than unity for the top (bottom) fibers, and 
increases slower (faster) than 0 and causes misalignment errors. As 0 
exceeds 0.8%, k/k approaches to unity and equation (29) leads to 

A0 ± j w A0 (30) 

The Implication of equation (30) Is that, at larger strain level, the relative 
error arisen from misalignment is always decreasing. 

The effect of strain rate can be further examined from Fig. 6, where k 
at the misdpan Is shown varying with the width of the specimen for two levels 
of 0. It can be concluded that the range of variation of lc at a cross-section 
Is wider and that a uniform distribution of strain rate Is achieved at a 
larger strain level for larger strain rate loading. 
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The misalignment error at the geometrical centerline (or center error) 
at the midspan is shown in Fig, 7 for three levels of strain rate. The 
center error, denoted by E (%) , is defined by fe - e)/e. It is seen that 

C C 

the center error is always decreasing in the i^astic range and vanishes at 
large strain. This result is also expected from equation (29). Generally, 
the center error is within 2% for annealed aluminum undergoing plastic defor- 
mation and vanishes mostly when the plastic strain is 2% or more. Although 
this investigation does not cover the effect of misalignment in the elastic 
range, it is expected that a result similar to that of Ref. [3] prevails, i.e., 
the misalignment effect is very significant near the interface between the 
elastic and plastic ranges, or the knee portion of the stress-strain curve. 

The error would then decrease as the plastic strain increases. 

The most important finding from Fig. 7 is that the misalignment effect 
reduces with increasing strain rate. An interpretation of this finding is 
that the material exhibits higher resistance to deformation at higher strain 
rate and this help retard the influence of misalignment. 

Misalignment errors of the top (E T ) and bottom (E^) fibers are plotted 
in Fig. 8. They have the same trend as the center error but with much greater 
magnitudes. The errors at lowest strain rate considered (k = 1) and at 0 = 0.1% 
are +50% and -45% for bottom and top fibers, respectively. Therefore, any 
investigation which is related to the local deformation of the extreme fibers 
at small plastic strain range may lead to invalid conclusion due to the signif- 
icant misalignment errors discussed above. 


Due to the limitation of computational funds, the curves for k = 0.75 and 
0.60 have not been completed at higher 9 level. However, it is expected 
that the curves will follow the trend shown by the dashed lines in Figs. 7-10. 
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The location of the N.A. Is plotted In Fig. 9 for the midspan. It Is 
seen that the smaller the strain rate is, the greater mobility the N.A. has, 
The peak values of d/h Is resided within the range of b equal to 0.5% to 1.0% 
for the three strain rates considered and can be approximately correlated to 
the points on the stress-strain curves where the stress-strain curves begin 
to flatten out. These points are shown by arrows In Fig. 2. 

It Is known from the discussion earlier that k/k ->• 1 when 0 > 0.8%. 

Under this condition equation (30) leads to 


0 ij " 0 


ij 



(31) 


where is a constant which varies along the cross-section. the defor- 
mation continues to the extent that the plas-xc strain dominates the total 
strain, the expression for the location of the N.A. can then be obtained from 
equation (7) and approximated by 



0 1B- 0 


(01 B - 0) - (0 1T - 0) 


(32) 


where the subscript 1 denotes the quantities at the midspan. Evidently, 
equations (31) and (32) suggest that the value of d^/h would reach a plateau 
(constant) as 0 is large. In fact, this is indeed the case in Fig. 9. For 
the case of k = 1.0, the d^/h curve reaches a plateau of value 0.51 for 0 
approximately equal to 1.75%. It is expected that the other two curves would 
also behave in a similar manner and find their own plateaus as 0 becomes large. 

The bending moment at the end of the specimen may be found from 

M = P(5* - La) 


(33) 
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This equation was derived In Ref, [3] and is shown in Pig. 3. The moment M is 

plotted in Pig. 10 for the three cases considered, It is found that the load 

1? and the distance d affect the end moment greatly. As far as the strain 

rate effect is concerned, the magnitude of the end moment is greater at greater 

strain rate. This is resulted from the larger stress that a specimen can 

sustain at a higher strain rate for a fixed plastic strain. Hence, the load 

P is also greater for greater strain rate. On the other hand, due to the 

influence of d , the end moment follows the trend of d and reaches a maximum 
s s 

value at the small strain range and then drops off rapidly and approaches a 
plateau at large strain. 

The deflection curves are shown in Fig. 11 for three levels of 0 for 
k ® 1 and two levels each for k = 0.75 and 0.60. Within the range of investi- 
gation, it is found that the rate of deflection decreases as the plastic strain 
Increases. 

Finally, the solutions presented above are compared with those obtained 
from a more restricted formulation given in Appendix A, in which the strain 
rate is assumed to be uniform throughout the specimen. Figures 12-14 show 
the results from both formulations for the case of k = 1.0. With respect to 
the solutions of the nonuniform strain rate formulation, the predicted mis- 
alignment errors and the shifting of the N.A. by the uniform strain rate for- 
mulation are overestimated at the small strain range and underestimated at 
the larger strain range. For the uniform strain r*ate formulation, the function 
k is uniform across the cross-section and the N.A. returns to the geometric 
centerline of the specimen at large strain. Although strain rate was not 
considered a factor in the work of Wu and Rununeler [3], it is seen that the 
formulation of [3] belongs to that of the uniform strain rate formulation and 
the N.A. would also return to the geometric centerline at larger strain level. 


The above two formulations are correlated during the course of numerical 
calculation. The results of the uniform strain rate formulation serve as the 
initial guess for the solution of the highly nonlinear equations obtained from 
the nonuniform strain rate formulation. The procedures of computation are 
presented in Appendix B. Since only the monotonic loading condition is of 
interest here, the step size AO is assigned for each iteration. It is found 
that for the purpose of rapid convergence, AG should increase with the plastic 
strain. However, due to the small number of elements used in the discrete 
model, the numerical results at larger strain level are unstable at the grip 
ends of the specimen (but is stable elsewhere in the specimen). It takes a 
great amount of trial and error procedure to select a correct step size so 
that the numerical output is convergent. It is anticipated that the problem 
of numerical instability mentioned above would be improved with greater number 
of elements taken in the discretized model. 
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5. CONCLUSIONS 

The following conclusions may be drawn from the present study of 1100-0 
aluminum specimens: 

(1) Based on the endochrcnlc constitutive equation, the formulation of 

the misalignment problem In the dynamic tension test has been developed. Two 
assumptions ore made: (a) plane cross-sections remain plane during deformation; 

and (b) the strain rate history effect is negligible in the computation of 
stress (equation (8)). The second assumption Is justifiable due to the 
extremely small variation In the strain rate throughout the tension test for 
each element of the discrete model, 

(2) The lower the strain rate is at the tension test, the more signifi- 
cant the misalignment errors become. Three levels of strain rate, l.e., 

1.30 x 10 '*» 1.24 x 10 and 7.63 x 10**^ sec ^ are investigated. In this 
strain race range, the misalignment error at the geometrical centerline of the 

specimen Is within 2% and vanishes mostly as the plastic strain Increases 

beyond 2 %. At the extreme fibers of the specimen, the misalignment errors 

approach to 50% at the very small plastic strain range. Therefore, any Inves- 

tigation related to the local deformation of the extreme fibers will have to 
account for the misalignment effect. 

(3) In the range of plastic deformation, the error introduced by mis- 
alignment decreases with the magnitude of plastic strain. Hence, it is con- 
jectured that the most significant effect of misalignment occurs at the inter- 
face between the elastic and the plastic range. This conjecture is in agreement 
with the finding of Wu and Rummler [3] that the most significant error occurs 

at the knee portion of the stress-strain curve. 


,*** ¥ 


(4) Ac che midspnn of Che specimen, Che variation of strain rate 
between extreme fibers decreases with increasing plastic strain. The strain 
rate reaches a uniform vulue at plastic strain equal to 0.8%. 

(5) The N.A. will shift toward the geometrical centerline of the 
specimen as the plastic strain increases. But, it will reach a limit and 
will not move bock to the geometrical centerline at large plastic strain as 
anticipated by the uniform strain rate formulation. 

(6) The misalignment errors and the shifting of the N.A. predicted by 
the uniform strain rate formulation in which the strain rate is constant for 
the whole specimen are always overestimated at the small plastic strain range 
and underestimated at the larger plastic strain range, 

(7) The formulation of Wu and Rummler [3] belong to the case of uniform 
strain rate, formulation. 
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FIGURE CAPTIONS 

Fig. 1 Strain Rate Sensitivity Function for 1100-0 At 

Fig, 2 Constant Strain Rate Stress-Strain Curves 

Fig. 3 (a) Specimen Dimensions and the Coordinate System, 

(b) The Load Train Configuration of the Symmetric Case of Misalignment 
Fig. 4 The Discrete Numerical Model 

Fig, 5 The Normalized Strain Rate Sensitivity Function of Top and Bottom 
Fibers at Midspan 

Fig, 6 Variation of k at the Midspan with Respect to the Width of the 
Specimen 

Fig, 7 Center Error at the Midspan 

Fig. 8 Misalignment Errors of the Tod and Bottom «■>->« 

Fig, 9 Location of N,A, at the Midspan 

Fig. 10 The End Moment 

Fig. 11 The Deflection Curves 

Fig. 12 Center Error for the Two Formulations (k - 1,0) 

Fig, 13 Errors at Bottom and Top Fibers for the Two Formulations (k * 1.0) 

Fig. 14 Locaton of N. A. for the Two Formulations (it * 1.0) 


APPENDIX A - Uniform Strain Roto Formulation 


In this formulation, the strain rate is assumed to be constant throughout 
the specimen In addition to the usual assumption that plane cross-sectlone 
remain pit as during deformation. The load and moment balance equations are 
obtained In closed form. Thus, In order to obtain a numerical solution, it 
la only necessary to eut the specimen into I element along the x-dlrectlon. 

Balance of force gives 



(Al) 


where subscripts T and B Indicate location at top and bottom fibers, respec- 
tively. Since, for monotonic loading, 


d5 « k(0)d0 



(A2) 


It follows 


de „ . 1 do 

dr, k(o) B dr a 


(A3) 


Note that in this formulation k(0) is constant throughout the test. 
Equations (1) , (Al) , and (A3) combine to give 


p m 

k 


j(f° + fa) 

2 (\26 2E/ 2 


+ (a 0 " °y> 



, X ( n x2 --2n+2 \ B 
+ 2E (o 0 " V 2 ) z T > 


for 5 > 0 


(M) 


z ** 1 + 65. 


where 


A2 


The balance of moment; le 

( 2 f C B 

M(x) - - I ay dA * -p eA + bp I ert dt (AS) 

A ; e T 

which may be written aa 

M(x) « - pe P + bp^CM^ + + M^) (A6) 

where 




A3 


and 


M, 


£ da 
r 2 dc 

h fp 


dC 


eV 


'0 -3 

n 2 


2/ 

o 0 (t, 0 


_ s"-n+3 , . .2 -*-2n+3 

V 2 + °0 (o 0 " V z 


(o n ~ 


V 


„ — 3n+3 

3 z 


(A10) 


T 


To arrive at the above equations, equation (A2) was integrated with the 
requirement of monotonic loading beginning at the zero stress state. Thus, 


e 



(All) 


Combining equations (All) and (1) , the following relation is obtained 


1 + Pke = 



1 

E 




(A12) 


Hence, the variables u, p, c^. . Eg, and Zg can be solved numerically from 
the system of equations (A4) , (A6), (A12) , (6), (22) and (25) subjected to the 
boundary conditions given by equations (27) and (28) . 


APPENDIX B - Numerical Procedures for the 

Nonuniform Strain Rate Formulation 


In this calculation, the specimen Is discretized Into 4 elements along 
both the longitudinal (I « 4) and the lateral (J * 4) directions. Algebraically, 
there are 55 unknown variables with 55 nonlinear equations. They are: 

Variables: o^ (20), e iT (5), e iB (5), (20), u i (4), d g (l) 

Equations: eq, (8) - 20, eq. (15) - 20, eq. (18) - 5, eq. (21) *-5, 
eq . (24) - 1, eq. (27) - 4. 


The computation of needs special attention. Specifically, the fol- 
lowing expression should be Incorporated into equation (10): 


An 


fti) - * ( V_y e ) - 

\ 0 / \ e - o/E / \ Ae - Ao/E / 


fE ( n ) _ E (n-1). (a M _ <n-l) w \ 

U i,j g ij y (q i,1 g l,1 j/C ( 

( -(n) _ i<nrl) ) _ Ca M „ 5 (n-l) )/E j 


(Bl) 


where Is the strain at the nth step and is that of the previous 

step. 

The numerical procedures are as follows: 

1. Use the solution of the uniform strain rate formulation as the initial 

guess to compute the first solution of the nonuniform formulation. This 

solution is the first loading step beyond the elastic limit. In obtaining 

Y Y 

this solution, the yield stress 0^ and the yield strain e^j are used as 

th 

the corresponding values for the (n - 1) step. In this case, equation 
(Bl) is reduced to 


B2 


(k 1 ) - *n HU- - -) (B2) 

\ 0 / V e - o/E / 

2. Increase 5 or 0 by a prescribed step size AC or A6. 

3. Calculate the mean stress, mean strain and load for the new step. 

4. Modify the true solution of the previous step according to the ratio 

and input the modified results as the initial guess for the 
current step of computation. 

5. Solve the system of equations and calculate the misalignment errors. 
Note that in this step, the expression (Bl) should be used. 

6. Repeat Steps 2 through 5. The computation terminates at a prescribed 
strain, Note that the step size A0 needs to be increased with 0. 


0s - 3.643 x 10 ' 2 
0r = 1.30 x io" 5 SEC' 1 
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